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Abstract 

Using the CLEO detector at the Cornell Electron-positron Storage Ring, we have measured the 
scaled momentum spectra, da/dxp, and the inclusive production cross sections of the charm mesons 
D~^, , D*~^ and D*^ in e"'"e~ annihilation at about 10.5 GeV center of mass energy, excluding 
the decay products of B mesons. The statistical accuracy and momentum resolution are superior 
to previous measurements at this energy. 
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I. INTRODUCTION 



We report the measurement of the momentum spectra of charged and neutral D and 
D* charm mesons produced at the Cornell Electron-positron Storage Ring, CESR, in non- 
resonant e~^e~ annihilation at about 10.5 GeV center of mass energy (CME) and observed 
with the CLEO detector. The and spectra each include both directly produced 
D's, and D's which are decay products of D excited states. From them we also derive the 
inclusive production cross section for these charm mesons. 

While very accurate data on bottom quark production from LEP and SLD have been 
pubhshed in recent years [l[ 0, 3L the data currently available for studies of charm 
fragmentation at 10.5 GeV CME p, |6|, are quite old and, by present standards, of poor 
statistical quality and momentum resolution. Our statistical sample is about 80 times larger 
than the our previous one and our current momentum resolution is about a factor of 2 
better. 

The spectra represent measurements of charm quark fragmentation distributions D^{x, s), 
i.e., the probability density that a c quark produces a charm hadron h carrying a fraction 
X of its momentum, ^/s being the "energy scale" of the process, the e^e~ CME in our 
case Experimental heavy-meson spectra in e^e~ collisions are important for theoretical 

and practical reasons: (i) they provide a component that is not yet calculable in predicting 
heavy flavor production in very high energy hadronic collisions, (ii) they can test advanced 
perturbative QCD (PQCD) methods, (iii) they can test the QCD evolution equations, and 
(iv) they provide information for best parametrization of the Monte Carlo simulations on 
which the analysis of many high energy experiments partially rely. 

Items (i) and (ii) are interconnected. The calculations of heavy flavor production cross 
sections in hadronic collisions [e.g., at the Tevatron and the LHC) are generally based on 
the factorization hypothesis, i.e., a convolution of (a) the parton distribution function for 
the colliding hadrons, (b) the perturbative calculation of the parton-parton cross section 
and (c) the parton fragmentation function D^{x, s). Items (b) and part of (c) (the parton- 
shower cascade) can be calculated, in the case of heavy quarks, using PQCD. Items (a) and 
the second phase of (c) (the hadronization phase) are intrinsically non-perturbative (long 
distance) processes: as of now, they must be provided by experiments. There is an ongoing 
theoretical effort to push the potential of PQCD to calculate the perturbative component 
of the fragmentation function. It needs tests and guidance from the experimental spectra of 
heavy flavored hadrons produced in e^e~ annihilation. De-convolving the calculated PQCD 
component from the experimental spectra, one obtains the non-perturbative component 
of the fragmentation function. Unphysical behavior of the result {e.g., negative values, 
extension beyond the kinematic limit) is indication that further refinement of the PQCD 
calculation is needed. Tests of this kind have been performed up to now on B production 
in e"'"e~ annihilation |^, |^ and in hadronic collisions 11, and on charm production 



in hadron and ep collisions Charm production in e"'"e~ annihilation provides 

a further testing ground of these theoretical attempts 3]. The larger value of Kqcd/t^c 
with respect to KQcn/^b makes these non-perturbative effects more evident than in bottom 
hadron production. 

Tests of the Altarelli-Parisi evolution equations 17, 3] have been performed by our 
collaboration with low sensitivity and over a relatively small energy interval, comparing 
the CLEO results with PETRA results. The spectra reported in the present paper can be 
compared with LEP results providing a test over the 10 to 200 GeV energy range. 
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Lacking rigorous calculations of the process of quark and gluon hadronization, QCD in 



spired Monte Carlo simulations have been built: the Lund String Model j20|, |21|, |2j| and 
Cluster Fragmentation j23||. These models have been implemented in Monte Carlo pro- 
grams (JETSET HI, UCLA HERWIG Q). In each case a number of parameters are 



introduced, to be determined by fitting the experimental distributions. Monte Carlo simu- 
lations of quark hadronization are used by experiments to determine detection efficiencies 
and to calculate some sources of backgrounds. The results presented here include a JETSET 
parametrization that produces spectra that agree quite well with the shapes of all spectra 
obtained in this analysis. 

In all these uses of our results, spectral shapes are most important, rather than the 
absolute cross-section values; therefore, shape is the main focus of our attention. 

In Sec.|n]we first list the charm mesons studied in our analysis along with the decay modes 
considered and then we describe the data sample analyzed and outline the procedures used to 
produce the spectra. In Sec. Illll we describe the Monte Carlo simulations we have generated 
and their use. In Sec. llVl we give details on how we extract the signal from the effective mass 
distributions, and in Sec.|V|we explain how the detection efficiency is estimated. Sec. I VII is 
devoted to discussing the checks we performed and the evaluation of errors. In Sec. IVIll the 
results, i.e., the charm meson spectra, are shown in the order given in Sec. |n] Our results 
for the inclusive production cross sections are given in Sec. IVIIII Our optimization of the 
JETSET parameters to reproduce our spectra is described in Sec. |TXl In two appendices we 
show plots of the detection efficiencies and provide detailed tables of the measured spectra. 



II. GENERAL ANALYSIS PROCEDURES 

We measure the momentum distributions of D^, D^, D*+and D*^ using the following 
decay modes (charge conjugates are implied throughout this paper): 

• D+ 



TT 



D* 



D*+ D%+ {K-n+)n+ 
D*+ D°n+ {K-n+TT+TT-)n+ 

We apply selection criteria to identify events with candidate D and/or D* that decay in 
one of these modes. We then extract the candidate D or D* mass distributions in twenty 
0.05 wide bins of the reduced momentum, Xp{D) = p/pmax, where Pmax (approximately 
4.95 GeV/c) is the maximum attainable momentum at the relevant beam energy. 
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We fit these mass distributions with appropriate signal and background functions. The 
distributions of signal yields vs Xp, corrected for detection efficiency, give the shape of the 
Xp spectra: the main goal of our analysis. We then divide these spectra by the integrated 
luminosity and the appropriate decay branching fractions to form the differential production 
cross section da/dxp for each channel. 

The use of different decay modes of the same meson provides a check on possible system- 
atic biases. 

The procedures used in the present analyses closely parallel those we used in measuring 
D and D* spectra from B decay. |26l | 

A. Data and Detector 

The e^e~ annihilation data sample used in this study was taken with the CLEO II. V 
detector 0, Q at CESR during 1995-1999. 

It consists of 2.9 fb~^ of the "continuum" (non-resonant) data sample at about 10.52 GeV 
CME (36 MeV below BBbar threshold) and the "0N4S" sample, comprising 6.0 fb"^ at 
10.58 GeV, the T(4S) peak. Assuming that the shape of the spectrum is the same at 
these two energies,^ we merge the two samples for charm mesons with momenta above the 
maximum kinematically allowed in B decay. For lower momenta we use only the continuum 
sample, thus reducing the statistics available in that region. All charm hadrons coming from 
B decays are thereby excluded. 

To combine the two parts of the spectra, Xp < 0.50 extracted from only the continuum 
sample, and Xp > 0.50 extracted from both the continuum and "0N4S" samples, using the 
well known 1/s dependence of the e~^e~ annihilation cross section into a pair of fermions 
(see Sec. 39 of ref. [3^, 2004), we scale the Xp < 0.50 spectra by the factor 



Here Cq and C4 are the integrated luminosities of, respectively, the "continuum" and "On4S" 
samples, and sq and S4 are the squares of the respective CMEs. The statistical sample for 
< 0.50 is a factor of three smaller than that for Xp > 0.50. 

The spectrum so obtained is then divided by the integrated luminosity, (£0 + '^4), and 
by the appropriate decay branching fraction to obtain da/dxp for each channel. 

B. Selection Criteria 

We select events using standard CLEO criteria designed to efficiently select e~^e~ an- 
nihilation into hadrons, while rejecting Bhabha scattering, e~^e~ —>■ fi~^fJ'~, and beam-gas 
interactions. At least three tracks are required. Events with three or four tracks must also 
have 65% of the center-of-mass energy deposited in the calorimeter. For those with five or 



Comparing our spectra with the corresponding ones at ^/s — 30.4 GeV |22j we estimated that the fractional 
difference between the D* spectrum at ^/s — 10.52 GeV and the one at 10.58 GeV is at most 0.075%, 
after normahzing one to the other. Because of this sample merging, our results effectively refer to CME 
Vi= 10.56 GeV 




6.0 (10.52)2 



(1) 
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more tracks the visible energy, summing both energy in tracks and neutral energy in the 
calorimeter, must exceed 20% of the center-of-mass energy. 

Tracks used to reconstruct a. D or D* are required to be the result of good tracking fits 
and to have an angle with respect to the beam line, 6, such that | cos ^| < 0.91. They are also 
required to be consistent with originating from the luminous region. Further, if they have 
momentum greater than 250 MeV/c, we require that the impact parameter with respect 
to the beam line be less than 3 mm, and that the distance between the point of closest 
approach to the beam line and the event vertex be less than 2.5 cm. 

We impose particle identification requirements based on specific ionization (dE/dx) and 
time of flight measurements for the track. The requirement is that the combined proba- 
bility of the chosen identification must be greater than 4%. 

Photon candidate showers detected in the central barrel region (| cos6'| < 0.707) of the 
crystal calorimeter are required to have a minimum energy of 30 MeV. Those detected in the 
forward calorimeters are required to have a minimum energy of 50 MeV. Photon candidates 
are also required to be well separated from the extrapolated position of all tracks, and 
the lateral shape of the energy distribution must be consistent with that expected from an 
electromagnetic shower. 

Candidate 7r° mesons are reconstructed from pairs of photon candidates. At least one 
of the two must be in the central barrel region. To improve the determination of the 7r° 
momentum, the two photon combination is kinematically fitted to the nominal mass. 
The combination is accepted if this fit has -P(x^) ^ 10%. The resulting 7r° 4-momentum is 
used in reconstruction. 



III. MONTE CARLO SIMULATIONS 

Monte Carlo simulations are used to estimate detection efficiencies. Continuum e~^e~ 
annihilation events are generated using the JETSET 7.3 [s^ package. The simulated events 
are then processed through a GEANT-based [lH simulation of the CLEO detector and 
reconstructed and analyzed as real data. 

The Monte Carlo simulations are also used for other purposes: (i) to provide a shape 
for the signal in the candidate D mass distribution fSec lIV|) . (ii) to estimate the D and D* 
momentum resolution fSec illl A]) , and (iii) to perform checks on the validity of our analysis 
procedures fSec HvHIVrai . 

We use two kinds of Monte Carlo simulations. In the first kind, the "signal Monte Carlo" , 
only e~^e~ —>■ cc events are generated at the JETSET stage, and an event is accepted only 
if the charm meson under study is present. That meson is made to decay only in the mode 
under study. The corresponding anti-charm hadron decays generically. We produce three 
signal Monte Carlo's, one for and two for for the two decay channels analyzed. The 
D's in these signal Monte Carlo's are the mix of directly produced D's and D's that are 
decay products of D*'s and other excited charm states. The mix is as generated by the 
physics simulation (JETSET). It follows that each one of these signal Monte Carlo's act also 
as signal Monte Carlo for D*'s decaying into that specific D channel. 

In the second kind of simulation, the "generic Monte Carlo" , all possible e"'"e~ hadronic 
annihilations are produced according to present knowledge 3^ • 



The three signal Monte Carlo's and the generic Monte Carlo accurately reproduce the 
D and D* signal shapes observed in data. Backgrounds in the signal Monte Carlo mass 
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distributions are much smaller than those in the generic Monte Carlo, which simulates more 
accurately the backgrounds in the data. 

Both kinds of Monte Carlo simulation are used to estimate the detection efficiency. For 
each D or D* meson and its decay chain, we find that the signal Monte Carlo and generic 
Monte Carlo-derived efficiencies are statistically compatible. This proves that the strong 
background reduction in the signal Monte Carlo does not affect the efficiency estimation or, 
vice versa, that the large background of the generic Monte Carlo introduces no appreciable 
bias in the detection efficiency. 

The two statistically independent Monte Carlo simulations allow internal checks of our 
procedures. We will refer to these as "generic Monte Carlo checks" . In a generic Monte Carlo 
check, we analyze the generic Monte Carlo as data, using the procedure to be checked. Then 
we correct the reconstructed momentum spectrum using the detection efficiency obtained 
from the signal Monte Carlo. Finally we compare this efficiency-corrected spectrum with 
the JETSET-generated spectrum that was the input to the generic Monte Carlo. This com- 
parison consists in calculating the of the bin-by-bin difference between the reconstructed 
and the input spectrum: 



where n is the number of bins, Ri and are the values of, respectively, the reconstructed 
and input spectra in bin i and 5Ri is the statistical error on Ri (the statistical errors on the 
input spectra are negligible). The resulting probability, or confidence level (CL), is the 
measure of the correctness of the analysis procedure being checked. If we normalize the two 
spectra to each other and recompute the the new CL is a measure of the correctness 
of our procedure in so far as the reconstruction of the shape of the spectrum is concerned, 
irrespective of normalization. 

In a generic Monte Carlo check, the comparison is with the input spectrum. It is sensitive 
to all sources of systematic error on the shape of the spectra, except for possible errors in 
physics and detector simulation, that are common to signal and generic Monte Carlo. Hence, 
insofar as the MC is correct, each check provides a comprehensive estimate of all systematic 
errors associated with the shape of the spectrum, for the procedure being checked. 

A. Momentum Resolution 

Comparison with theoretical calculation may involve the moments of the spectra: 
J^x^^dx. In order to minimize correlations between adjacent Xp bins, the Xp bin size 
should be chosen to be substantially larger than the Xp resolution. It is then important to 
know the momentum, and hence the Xp, resolution in our analysis. Using the CLEOG Monte 
Carlo simulation which reproduces rather accurately our track and shower measurement 
errors, we plot the difference between the reconstructed Xp and input Xp (from JETSET). 
Fig. ^ shows this resolution distribution for the mode D^^K^tt^ for all momenta. The full 
width at half maximum (FWHM) is 0.008, i.e., 16% of the bin size (0.050). The resolution 
(FWHM) varies monotonically with momentum, from 4% of bin size at Xp = 0.10 to 18% 
for Xp = 0.95. For the other channels the resolution is likewise a small fraction of the bin 
size. 




(2) 
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FIG. 1: Resolution in Xp for the D^^K~it~^ channel. All momenta. 

IV. CANDIDATE MASS DISTRIBUTION FITTING 

For the and analyses we select candidate daughters, add their four-momenta, and 
calculate the invariant mass Mcand of the charm meson. Multiple candidates in the same 
event are accepted. 

In the D* case we obtain the Mcand distribution for the D*^ associated with the D* by 
selecting D* candidates with Q = M*^„^ — Mcand — f^-n in the signal region for D* decay. 
Here M*^„^ is the invariant mass of the decay products of the candidate D* . Random D-tt 
associations are subtracted using the Mcand distribution for events in the side bands of the 
D* signal in the Q distribution.^ 

Fig. El shows examples of the Mcand distributions for three different D* decay modes, 
for events with Q in the signal region and for those in the Q side bands. The residual 
background after the subtraction is due to D candidates from random track association. 



^ The signal and the side-band regions are defined as follows. We fit the "global" {i.e. all momenta) Q 
distribution with a Double-Gaussian plus suitable background. The ratio SIG2/SIG1 of the widths of the 
two Gaussians is, in all cases, about 2.2. We choose the signal region to be MEAN±n*SIG2, where n 
(that turns out to be about 2 in all channels) is evaluated from the Gaussian Integral tables, requiring 
that the whole area of the narrow Gaussian plus the area within ±n*SIG2 of the wider Gaussian result 
in a 98% of the Double-Gaussian area. For the side bands, on each side, we skip n*SIG2 and then take a 
region n*SIG2 wide. 
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FIG. 2: Examples of M^and distribution for two D*'^ decay channels and one of the D*^ channels 
analyzed, (a) D*+ ^{K~-k+)t:+ . (b) L>*+^(K-7r+7r+7r-)7r+. (c) L»*o^(K-7r+7r+7r-)7rO. They 
show the Mcand distribution for Q in the D*^ signal region and for Q in the side bands. 



The choice of the signal shape used to fit the M^and distribution was studied and discussed 
in detail in a previous paper j2a]. A Gaussian function does not give a sufficiently accurate 
parametrization of the D signal. Track measurement errors vary because of the geometrical 
orientation of the D decay products in the detector, because of different momenta of the 
decay tracks and overlap with other tracks. That study concluded that a satisfactory choice 
for the D signal shape is a double Gaussian, i.e., the sum of two Gaussians constrained 
to have the same mean. A different choice of signal fitting function is the signal shape 
obtained from the Monte Carlo simulation where, for each track, we can identify the input 
particle that generated it. We call the signal mass histograms thus obtained (one for each 
momentum bin) the "TAGMC shape" . To compare these two choices we repeat a test that 
was performed in the previous paper ji^, on the D^^K~t[^ channel, as follows. 

We repeat the data analysis, replacing the double Gaussian with the TAGMC shape. 
With this signal shape we obtain excellent fits, although not superior to the double Gaussian 
fits. We use MINUIT to find the compatibility of the two spectra. We fit one using the 
other as fitting function. The fitted relative normalization parameter is 1.016 ± 0.007, and 
the CL of the fit is 93.8%. The two spectra are compared in Fig. |2fa) after normalizing one 
to the other. To find if there is any Xp dependence of the difference between the spectra 
obtained by the two methods, we took the bin-by-bin fractional difference between the two 
spectra (Fig. IHl^b)) and fitted it to a constant, resulting into a CL=91.0%, consistent with 
no difference between the two choices of signal shape. The results obtained using the double 
Gaussian as signal shape, are compared with the TAGMC shape to estimate the systematic 
error on the total cross sections due to the uncertainty on the signal shape. 

The suitability of the double Gaussian as fitting function is also confirmed by the goodness 
of the fits: in all the channels, the fit confidence levels are evenly distributed between 0.0 and 
1.0, as they should be. A quadratic polynomial is used to fit the combinatoric background 
in each of the seven channels. 

The fits of the Mcand distributions are over the whole 1.70-2.02 GeV range shown in the 
figures, except for the K^ii^n^ case, where we exclude the 1.96-2.02 GeV (D*^) 

region, and for the K^tt^ case, as explained in the next subsection. The fitted area of 

the double Gaussian (or the result of the COUNT procedure described in Sec. IIV 21 below) 
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FIG. 3: (a) Overlay of spectra (data) from double Gaussian and TAGMC shape signal fitting; 
(b) fractional difference of the two spectra. 



is the "raw" yield for that Xp bin. 

In the next two subsections, we discuss additional backgrounds in the Mcand distribution 
from the K^tt^ channel, and describe an alternative procedure, the COUNT method, 

to estimate the raw yield in the K' 'TT+TT 7r+ channel. 



1. The D° K-TT+ case 

In the — > K^n^ case (direct or from D* decay) additional backgrounds must be 
considered: D'^ decays to K~K~^, vr^Tr^, K~ and ^ ^K~^7r~ misinterpreted as K~7t~^. 
The shapes of their Mcand distributions are obtained from Monte Carlo simulation. 

The K^p^ background is very small and contributes only to the 1.70 < M{K^7i^) < 
1.75 GeV mass region. This contribution is excluded by not considering this mass region in 
the fit. 

The background due to Kit switched identities shows as a very broad enhancement cen- 
tered at the signal position. For Xp > 0.20, this enhancement is so broad that it can be 
easily accommodated by the quadratic term of the polynomial background function. For 
small Xp, it is narrower, but contributes negligibly. The amount of this background is fixed 
to a momentum dependent fraction determined by Monte Carlo simulation. 

The backgrounds due to decays to K~K^, tt^tt^ do not contribute to the peak, 
but, if ignored, would result in a very poor fit of the background. Such a fit overestimates 
the amount of background under the signal and thus underestimates the amount of signal. 
The D^^K~ background level is a parameter to be fitted. Because of lack of statis- 
tics, the amount of D^^ti^tx^ background is constrained to a fixed fraction (0.357) of the 
D^^K^ background, based on the known relative branching ratio [s^. The tttt contri- 
bution is very small, and alternative methods of accounting for it cause negligible changes 
in signal yields. 

Fig. m shows data in three representative momentum intervals, demonstrating how the 
background is built up from the four contributions. All four background components are 
needed to extract the yield. 
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FIG. 4: Buildup of the background from its components to fit the M{K vr^) distribution. The 
sohd histogram is data. Notice the offset on the yield axis. 

2. The K-TT+TT+TT- case: the COUNT method 



In the case of the K^ti^ti^ii^ decay, direct or from decay, in addition to using 

a double Gaussian as fitting function for the signal, we use a different procedure that leads 
to results that are statistically competitive. In the — > K~n~^7r~^7r~ case, the signal is 
quite narrow and the background is smooth over a wide M^and region. We exclude the signal 
region and fit the background to a polynomial. The signal region is centered on the mean of 
the double Gaussian fit and its range is chosen so as to contain the entire signal. We then 
count all events in the signal region and subtract the background obtained from this fit. The 
result of this subtraction is the measured signal yield. We perform this procedure on data for 
three choices of the signal region: 1.810-1.920 GeV, 1.820-1.910 GeV and 1.830-1.900 GeV. 

We repeat this procedure on the generic Monte Carlo, thus performing the generic Monte 
Carlo check, described in Sec. HTH The 1.820-1.910 GeV exclusion gives the best CL: 28%. 
The narrower exclusion gives the worst CL: 6%. The wider exclusion gives an acceptable 
CL: 22%, in part, because the wider the exclusion region is, the larger the statistical error 
becomes. Based on these results, we choose the data spectrum obtained with the 1.820- 
1.910 GeV exclusion as our result. The bin-by-bin rms spread of the three data spectra 
obtained with different signal region exclusions is taken as the estimate of the systematic 
error of this procedure. 

We have two valid measurements, one from the COUNT method and the other from 
double Gaussian fitting of the signal, both performed on the same statistical sample. Hence 
we take as result the bin-by-bin arithmetic average of the spectrum obtained by double 
Gaussian fitting and the one obtained by the COUNT method with the optimal choice of 
the signal region exclusion: 1.820 < M^and < 1-910 GeV. 



A. Fit parameter smoothing 

The shape parameters of the signal and background functions are expected to depend 
smoothly on Xp. By imposing this smoothness of the shape parameters we suppress, in part, 
the bin-to-bin (in Xp) statistical fluctuations in the spectra. This improves the accuracy 
of the shape of the spectra, particularly at low Xp where statistics are poor. This param- 
eter smoothing procedure was used also in our measurement of charm meson momentum 
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spectra from B decay |26|. In the last paragraph of this subsection we show the extent of 
improvement obtained. 

The parameters considered are: the mean of the double Gaussian (common to the two 
Gaussians), the width of the narrower Gaussian, cxi, the ratio of the widths of the wider to 
the narrower Gaussian, o"2/o"i, and the ratio of the area of the wider Gaussian to the total 
area, A2/Atot- We impose this smooth behavior by fitting the Xp dependence of each shape 
parameter to a polynomial, at most quadratic, in Xp. 

We proceed in stages. We start by smoothing the parameter that shows the least fluctu- 
ations and repeat the Mcand distribution fitting for all the Xp bins, fixing that parameter to 
the value given by the smoothing function. We do this in sequence for all shape parameters. 
If a parameter does not show appreciable statistical fluctuations, we may skip smoothing it. 
It may take up to five iterations to smooth all the parameters. 

At each stage we get a new Xp spectrum and check that we have not introduced any 
distortion to that spectrum. The check is performed by calculating the bin-by-bin ratio of 
the new spectrum to the original one where all the parameters were allowed to float (the "no 
smoothing" spectrum). This ratio should show only random fluctuations around unity. If 
the ratio shows any trend vs Xp, e.g., if a slope and/or a curvature is needed to describe the 
Xp dependence of the ratio, that smoothing stage is discarded. Fig. El shows three examples 
of these checks. When we perform a fit of the ratios to a constant function (=1), we 
obtain CL of 94.6%, 91.0% and 38.0% respectively for the three examples shown. These are 
typical for all the retained smoothing steps. 
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FIG. 5: Ratios of data spectrum after double Gaussian shape parameter smoothing to the one 
obtained without smoothing: (a) D+^K-Tr+TT+, (b) D°^K-tt+, (c) D*+^{K-TT+)Tr+ . 



We perform the smoothing procedure varying the sequence of smoothing stages. Each 
change of sequence leads to a spectrum that is slightly different from the other ones. If the 
CL of the generic Monte Carlo check for one of the sequences is considerably higher than 
the CL for the other ones, we take that spectrum as our result. 

Comparison of spectra derived from different smoothing sequences provides a measure of 
the associated systematic error, as explained in Sec. I VI CA 

We use the generic Monte Carlo check discussed in Sec. Illll to see if the smoothing pro- 
cedure improves the agreement between the reconstructed and the original spectrum, i.e., 
the spectrum that is the input to the Monte Carlo simulation. In the D^^K~tt^ case, 
when there is no smoothing, the spectrum produced by the analysis fits the original ("true") 



12 



spectrum with a = 25.1 for 15 d.o.f., i.e., CL = 5%.^ When smoothing is used, the 
spectrum produced by the analysis fits the original spectrum with = 7.0 for 15 d.o.f., 
i.e., CL = 95%. Thus, in this case, parameter smoothing produces a dramatic improve- 
ment. In the case of D^^K'tt'^tt'^ , the CL improves appreciably from 7% to 13%. In the 
D*^^{K^ 71^)71^ case, where the CL is already 93% without parameter smoothing, there 
is only a small improvement to a CL=97%. In the D*^^{K~ n~^)7r^ case the improvement 
is from CL=59% to CL=75%. As expected, the improvement is strong when the initial set 
of parameters show large fluctuations, smaller when the parameters show a fairly smooth 
behavior to start with. 

V. DETECTION EFFICIENCY 

For each channel we have two independent and statistically-compatible estimates of the 
detection efficiency, as explained in Sec. lIIII We take their weighted average, thus appreciably 
reducing the statistical error on the detection efficiency. 

The detection efficiency should be a smooth function of Xp. We use a second order 
polynomial to fit the Xp dependence of the detection efficiency averaged over the signal and 
generic Monte Carlo. Adding a cubic term does not improve any of the fits. We call the 
result of this fit the "smoothed efficiency". In Appendix^ we show the detection efficiency 
dependence on Xp for all the mesons and decay modes analyzed. In Figs. |22l 1211 12S1 the 
detection efficiencies obtained from the signal and generic Monte Carlo's are plotted, and 
the curve resulting from the fit of their average to a polynomial is overlaid. This procedure 
results in a strong reduction of the statistical errors on the detection efficiency. 

The detection efficiency corrected spectrum is obtained by dividing the raw signal yield 
by the smoothed efficiency, bin-by-bin in Xp. 

VI. CHECKS AND ERROR ESTIMATION 
A. Two Checks 

1. Generic Monte Carlo checks 

For each procedure used to reconstruct the spectra, we perform a "generic Monte Carlo 
check" , as described in Sec. Illll The confidence levels reported below in Table HJ show the 
consistency of the reconstructed spectrum with the original one. Since our interest is in the 
consistency of the shapes of the two spectra, we do the comparison after normalizing the 
areas of the the two spectra to each other. The normalization differs from unity by at most 
2.6%. Notice that in the generic Monte Carlo checks we can only use the signal Monte Carlo 
efficiency, not the averaged, smoothed efficiency described in the previous section fSec. IV|). 



^ Since our aim is to measure the shape of the spectra, irrespective of normalization, these 8.nd related 
CLs are calculated after normalizing the reconstructed spectrum to the original one, thus resulting in an 
increase of the CLs. 
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TABLE I: Confidence levels of the fit of the generic Monte Carlo reconstructed spectrum to its 
input spectrum for the seven decay channels analyzed. 



Decay channel 


C.L. 


Decay channel 


C.L. 


Decay channel 


C.L. 


D+ K-Tr+-K+ 


18% 




72% 


K-TT+TT+TT- 


56% 


D*+ {K-7r+)TT+ 


70% 


D*+ (K-7r+7r+7r-)7r+ 


37% 




76% 






Z)*o ^ (^-^+^+^-)^o 


99% 







2. Comparison of spectra from different decay modes 

In the D^, D*~^ and D*^ cases we obtain the respective spectra from two different decay 
modes. We checked that the spectra from the two different decay modes are statistically 
compatible. We calculate the of difference, using only the statistical errors. The 
corresponding confidence levels are, respectively, 28%, 100% and 0.09%. After normalizing 
one to the other the confidence level become: 85%, 100% and 84%. This test, however, is 
not very stringent because the comparison is dominated by the large statistical errors of the 
D'^ K~7:^t:^t:~ channel. 

B. Statistical Errors 

The statistical errors on the efficiency-corrected yields are obtained by adding in quadra- 
ture the statistical error on the raw yield and the statistical error on the smoothed efficiency 
(Sec. E}. The latter is generally considerably smaller than the former. 

C. Systematic Errors 

We discuss here systematic errors that could affect the shape of the differential cross 
section da/dxp, although some of them are found to be independent of Xp. Additional 
systematic errors that affect the normalization of the differential cross section, but not its 
shape, will be discussed in Sec. IVIIll on total cross sections. 

1. Errors found to be independent of Xp or negligible. 

We consider the following possible sources of systematic errors: (1) the choice of signal 
fitting function, (2) possibly incorrect simulation of the initial state radiation, (3) effects of 
swapping between background curvature and width of the wide Gaussian in M^and distribu- 
tion fitting, and (4) effects of low detection efficiency for very low momentum tracks. 

The test, described in Sec. IIVI that uses a signal fitting function other than a double 
Gaussian, gives us a measure of the sensitivity of our results to the choice of signal fitting 
function. Based on that test, we attribute a systematic error of 1.6% from the choice of 
signal fitting function. The test shows no momentum dependence of the difference between 
the two methods. 
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We have considered the possibihty that inaccurate simulation of initial state radiation 
(ISR) may have introduced a systematic error in our estimate of the detection efficiencies. 
We compare the detection efficiencies discussed in Sec. |3with those obtained from Monte 
Carlo events where no ISR was produced. As expected, the latter is slightly higher than 
the former, but only by 1.1%, and its dependence on Xp is negligible. Since our Monte 
Carlo does simulate the initial state radiation, the uncertainty is only in the accuracy of the 
simulation. We thus take half of that, 0.5%, as contribution to the systematic error on the 
cross sections. 

Since the momentum dependence of these two uncertainties is found to be negligible, we 
take them into account only as errors in the total cross sections (Sec. IVIII|) . 

We considered the possibility of swapping between a background that is highly curved in 
the signal region, and the wide component of the double Gaussian. The only two channels 
that show an appreciable background curvature are D''^i^'^7r+7r+7r^and D^^K^-k^-k^ . In 
the first case the full compatibility of the fits with the results of the COUNT procedure 
(subsect. II V 21 CL> 96% for both Monte Carlo's and for data), shows that this swapping, 
if it exists, generates an error much smaller than the statistical error. In the case we 
performed the same test with the same result. 

We considered the possibility of errors in the D*'^ detection efficiency because of the 
very rapid decrease in the charged track detection efficiency for momenta below 120 MeV/c. 
The detection efficiency is practically zero below 70 MeV/c.^ We studied in detail the 
momentum distribution of the charged vr^ daughter of the D*^ (the "slow pion") as a 
function of Xp(D*^). Only for Xp{D*^) < 0.40 are there slow pions with momentum below 
120 MeV/c. From the momentum dependence of the track detection efficiency and the D*^ 
isotropic decay distribution IsS^], we can calculate the D*"^ detection efficiency. The result is 
consistent with the one resulting from our generic and signal Monte Carlo simulation within 
their statistical errors. 

Since we find the errors from these last two sources to be negligible, we disregard them. 

2. Errors that affect the spectra shapes 

The different sequences of parameter smoothing stages (described in Sec. IIV A|) lead to 
slightly different resulting spectra. We calculate the root-mean-square (rms) spreads of the 
yields for each Xp bin over the spectra from different sequences. Since these rms spreads 
fluctuate statistically from bin to bin, as expected, we average them over groups of three 
bins. We take these rms spreads as systematic errors on the yields. 

As stated in Sec. lIIIl we have both generic and signal Monte Carlo samples of events, and 
to the extent that our Monte Carlo correctly simulates data and detector, we can perform 
a test which give comprehensive information on all systematic errors associated with our 
analysis procedures. We take the bin-by-bin difference between the generic Monte Carlo 
reconstructed spectrum and the input spectrum, and divide this, bin-by-bin, by the input 
spectrum, resulting in the distribution of the fractional difference vs Xp. The weighted av- 
erage, over the entire Xp range, of the absolute values of these fractional differences (where 
the weights are the inverse square errors on the differences) can be considered as an esti- 



^ The charged track detection efficiency has been carefuUy studied in a series of CLEO internal documents 
(unpubhshed) . 
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mate of the systematic error. It varies from 0.6% for the D^—^K~ti^ channel to 1.4% for 
the D^^K~TT^Ti^ channel. The distributions of the fractional differences show negligible 
dependence on Xp, meaning that this estimated systematic error does not seem to affect 
the shape of the spectra. Nevertheless we include these average differences as a component 
of the systematic error on the measured yields. In principle, this estimate of the system- 
atic error takes into account also the "rms spreads" discussed in the previous paragraph. 
We decided, however, to be conservative, and have combined them in quadrature to obtain 
the total systematic error. Even with possible overestimate, generally the systematic error 
makes the total error larger than the statistical error by only 10% to 30%. 

D. Total errors 

The statistical errors and the two systematic errors affecting the spectra shapes are listed, 
channel by channel, in Table |3 - IXII in Appendix El These three errors are combined in 
quadrature to give total errors relevant to the shape of our spectra. 

VII. RESULTS ON THE SHAPE OF THE SPECTRA 
A. The Final or Combined Spectrum. 

For each D or D* meson and its decay chain, we obtain the spectrum fitting the signal 
with a double Gaussian after smoothing the Xp dependence of the Gaussian parameters, 
as described in Sec. IIV Al When we also employ the COUNT method, as explained in 
Sec. IIV 21 the spectrum that we report is the average of the spectrum obtained by fitting 
a double Gaussian and that obtained with the COUNT method. Details specific to each 
channel, are given in the sections showing the respective spectra. 

The spectra shown in the following are differential, inclusive pro- 
duction cross sections, da{e'^e~ D^*^X)/dxp at y/s = 10.58 GeV 
fully corrected for detection efficiency and decay branching ratios. We 
use the following decay branching ratios: B{D^^K~7t^) = {3.82±0.09)%, 
B{D'^^K-n+7r-n+) = {7A9±0.31)%, B{D+^K- 11+ n+) = {9.0±0. 6)%, 
B{D*+^D°7r+) = {Q7.Q±0.5)%, B{D*^^D^tt°) = {Q1.9±2.9)%. They affect only the 
normalization, not the shape, of the spectra. Uncertainties in the branching ratios will be 
reflected in the systematic errors on the total cross sections. Sec. IVlllI 
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B. Spectrum 

Fig. ini shows examples of fits to the Mcand distributions in three representative Xp bins, 
using fully smoothed parameters. Our result is shown in Fig. |71 and tabulated in App. IHl 
Table |V1 The spectrum shown is obtained after smoothing the Xp dependence of the double 
Gaussian shape parameters (see Sec. IIV A|) using the sequence that gives the best CL in the 
generic Monte Carlo check fSec. IVlX|) . 
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FIG. 6: Three examples of M{K tt^t:^) distribution fits. Notice the large vertical scale offsets. 
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FIG. 7: Differential cross section da{e^e~ — > X) / dxp in pb from the D'^^K^tt^ti^ decay 
mode, (a) shows explicitly the total and statistical errors, (b) the same spectrum overlaid with 
the JETSET spectra generated with two different sets of parameters (Sec. lfX() . 
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C. D'^ Spectrum 

1. Spectrum from D^^K^n^ 

Fig. ISl shows examples of fits to the Mcand distributions in three representative Xp bins, 
using fully smoothed parameters. 
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FIG. 8: Three examples of M{K^tt^) distribution fits. Notice the large y offsets. 



The inclusive, differential production cross section obtained from this decay mode is 
shown in Fig. Eland in App. ^ Table IVTl It is obtained after smoothing the Xp dependence 
of the double Gaussian shape parameters (see Sec. IIV A|) using the sequence that gives the 
best CL in the generic Monte Carlo check (Sec. IVL^ . 




0.25 0.50„ 0.75 1.00 



FIG. 9: Differential cross section da{e'^e — > D'^ X^/ dxp in pb from the — > K ■k'^ decay mode. 



2. Spectrum from D^^K -k'^tt'^t: 

Fig. Uni shows examples of fits to the Mcand distributions in three representative Xp bins, 
with no parameter smoothing. Because of the large statistical errors, we find the Gaussian 
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parameter smoothing procedure to be umeliable. However, as discussed in Sec. IIV 21 for 
this mode we use also the COUNT method with three different widths of the excluded signal 
region in order to get part of the systematic error on this procedure. 



D°^K Jt+Jt"3T+ 



61000 



51000 



29000 



1700 



*1 630903-008 




1 1 M 1 M M 1 M M 1 1 

:0.90< X <0.95 


' 







1.70 1.80 1.90 2.00 1.70 1.80 1.90 2.00 1.70 1.80 1.90 2.00 

M (K"jt+3t"jr+) (GeV) 

FIG. 10: Three examples of M(K^Tr^TT^Tr^) distribution fits. Notice the large y offsets. 
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FIG. 11: Differential cross section da{e~^e D^X)/dxp in pb from the K Tr+Tr+vr decay 

mode. 



The inclusive, differential production cross section obtained from this decay mode is 
shown in Fig. ^2 and tabulated in App. ^ Table IVlll It is the arithmetic average of the 
one obtained by double Gaussian fits (without any Gaussian parameter smoothing) and the 
one produced with the COUNT procedure, excluding from the background fit the 1.820- 
1.910 GeV region. For the final statistical errors we take the average of the statistical errors 
associated with the two methods. 
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3. The Average Spectrum 



The weighted average of the spectra obtained from the two decay modes analyzed is 
shown in Fig. IT^ and tabulated in App. IHl Table IXITl The two JETSET generated spectra 
are explained in Sec. IIXI 
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FIG. 12: Differential cross section da{e'^e~ — > D^X)/dxp, weighted average of the spectra from 
the D^^K~7r~^ and D^^K~7r~^7r~^7r~ decay modes, overlaid with the JETSET spectra generated 
with two different sets of parameters fSec. IIX|) . 



D. The D*+ Spectrum 



In Sec. |TV| we described our procedure for selecting D*^ candidates. The difference 
between the two Mcand distributions shown in Fig. |21 eliminates random D^n^ associations. 

1. D*+ Spectrum from D*+^D°TT+^(K-7r+ >+ 

The subtracted Mcand distribution (Fig. ITSj) shows the additional backgrounds present in 
this decay mode. They have been handled as described in Sec. IIV 11 
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FIG. 13: Three examples of fits of the subtracted M{K^tt'^) distributions for D* 
{K~7r~^)7r^ candidates. 
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The spectrum is shown in Fig. and tabulated in App. El Table I Villi It is the one 
obtained after smoothing the Xp dependence of the double Gaussian shape parameters (see 
Sec. IIV A|) using the sequence that gave the best CL in the generic MC check (Sec. IVI A|) . 
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FIG. 14: da{e+e- D*+X)/dxp, from the D*+ D^ti+ {K-tt+)ti+ decay mode. 



2. D*+ Spectrum from D*+^D^ti+^(K-tt+tt+tt- )-k+ 

Just as in the case of K~7t~^tt~^7t~ , taking advantage of the narrowness of the 

signal over a background that is smooth and well determined over a large region, we use the 
COUNT procedure described in Sec. IIV 21 with the signal region exclusion as optimized in 
that analysis (1.820-1.910 GeV). The Q selection reduces drastically the background with 
respect the case, and we obtain good double Gaussian fits of the signal as shown, for 
three representative Xp bins, in Fig. IT^ 
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FIG. 15: Three examples of fits of the subtracted M^K^ir^ir^ir^) distributions for D* 
D^TT^ {K^TT^TT^TT^)'ir^ Candidates. 
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FIG. 16: d(j{e+e- D*+X)/dxp from the D*+ 0^71+ {R- 7:+ 7^-71^)11+ decay mode. 



The spectrum is shown in Fig.^Jand tabulated in App.El Table HXl It is the arithmetic 
average of the one obtained by double Gaussian fit, after full smoothing of the Xp dependence 
of the double Gaussian shape parameters (see Sec. IIV A|) . and the one produced with the 
COUNT procedure, excluding from the background fit the 1.820-1.910 GeV region. 



3. The Average D*~^ Spectrum 
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FIG. 17: Differential cross section da{e~^e~ D*^X)/dxp, weighted average of 

L>*+^D°7r+^(K-7r+)7r+and D*+^L>°7r+^(/r-7r+7r+7r~)7r+spectra, overlaid with the JETSET 
spectra generated with two sets of parameters fSec. IIX|). 

The weighted average of the spectra obtained from the two decay modes analyzed is 
shown in Fig. [T7I and tabulated in App. iBl Table IXlTl The two JETSET generated spectra 
are explained in Sec. IIXI 
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E. D*" Spectrum 



To suppress random D^tt^ associations, we use the subtraction procedure already used 
for the D*~^ cases and illustrated in Fig. |21 
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FIG. 18: Three examples of fits of the M{K vr^) distributions for D*^ 
candidates. 
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Fig. ITHl shows three examples of fits of the subtracted M^and distribution for this channel. 
Here too we add to the fitting functions the backgrounds described in Sec. IIV II 

The differential cross section is shown in Fig. and tabulated in in App. El Table IXl 
Among the different stage sequences in smoothing the Gaussian parameters (see Sec. IIV A|) 
we choose the one that gives the best CL in the generic MC check (Sec. IVI A|) . 
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FIG. 19: do-(e+e- ^ D*^X)/dxp, from the D*^ Z)°7r° ^ {K-TT+)Tr^ decay mode. 
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2. D*^ Spectrum from D*^-^D^tt^-^(K-'k+tt+'k- )tt' 



Fig. HOI shows, for three representative Xp bins, the fits of the subtracted Mcand distribu- 
tion, using a double Gaussian and a polynomial background. 
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FIG. 20: Three examples of fits of the subtracted M{K tt^tt tt^) distributions for D*^ D^tt^ 
{K~7r~^'K~7:'^)TT^ candidates. 



Because of the smaller decay branching ratio and the smaller detection efficiency, due to 
the presence of a 7r°, the statistical errors are quite large, especially for Xp < 0.50, where 
we can use only the continuum events. We have used both the COUNT procedure and the 
double Gaussian signal fitting (without parameter smoothing) to get the yield. 
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FIG. 21: da{e+e- D*^X)/dxp, from the D*" D^tt" ^ (i^-7r+7r-7r+)7r" decay mode. 

The spectrum is shown in Fig. |^ and tabulated in App. El Table IXII It is the arith- 
metic average of that obtained by fitting the signal with the double Gaussian (smoothed 
parameters) and the one obtained by the COUNT method using the signal region exclusion 
optimized in that analysis (1.820-1.910 GeV). 
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3. The Average D*^ Spectrum 



The weighted average of the spectra obtained from the two decay modes analyzed is 
shown in Fig. |221and hsted in App. IHl Table IXITl The two JETSET generated spectra are 
explained in Sec. IIXI 
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FIG. 22: da{e+e- D*^X)/dxp, weighted average of the D*o^D°7r°^(i^-7r+)7r° and 
D*^^D^-K^^[K~-K^iT'^T:~)T:^decay modes. Overlaid are the JETSET spectra generated with two 
sets of parameters (Sec. IIX() . 



VIII. RESULTS FOR THE TOTAL CROSS SECTIONS AND AVERAGE Xp 

The production cross section for each channel is shown in Table lllll It is calculated 
by summing each differential cross section bin-by-bin. The first error in the table is the 
statistical error, obtained by combining in quadrature the statistical errors in each bin. If 
the yield in the lowest few bins cannot be reliably measured, the cross section is corrected by 
extrapolating the spectrum to Xp = using the JETSET distribution that fits the spectrum, 
discussed in Sec. IIXI This correction is between 0.2% and 6%. 

In Table HTl we list, channel by channel, the components of the systematic error on the 
production cross sections. In the first column we report the rms spread of the cross sections 
obtained by the four or five smoothing sequences used for each channel. The discrepancy 
between the areas of the input and reconstructed spectra in the generic Monte Carlo check 
fSec. lVlX|) . is shown in the second column. In the third column we list the percent difference 
between the integral of the spectra obtained using the double Gaussian and the one that uses 
the TAGMC signal shape fSec. IIV|) . This error is not considered for the channels where the 

decays to /T'Tr+Tr+Tr", because of the use of the COUNT procedure for those channels. 
We assume a 10% error on the extrapolation and show it in column 4. The remaining 
systematic errors are estimated and discussed in a series of CLEO internal notes and are used 
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in all CLEO analyses where they are relevant. We estimate a 1% per track uncertainty in 
the charged-track detection efficiency and 0.8% per track for particle identification efficiency. 
The choice of track quality and geometrical cuts result in an error of 0.5% also per track. 
The per track errors, being coherent, are multiplied by the number of tracks in the decay, 
and are shown in columns 5, 6, and 7. The vr" detection uncertainty is estimated to be 3% 
per 7r° (column 8). As discussed in Sec. I VI CI we attribute a 0.5% error due to possible 
inaccuracies in the Monte Carlo simulation of the initial state radiation. The error on the 
integrated luminosity is estimated as 1.9%. 



TABLE II: Systematic errors described in the text. Some are listed as percent of the cross section, 
other ones directly in pb. The momentum dependent systematic errors are listed also in the tables 
in App.^ The error due to the uncertainty on the branching ratio is shown only in Table HTll 
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These systematic errors are combined in quadrature to give the systematic error on the 
cross section, the second entry in Table UTTl 

TABLE III: Total production cross sections and average Xp, as derived from each decay mode. The 
cross section errors are, in this order, the statistical error, the systematic error and the error due 
to the uncertainty on the branching ratio. 



Decay channel 


Total Cross Section (pb) at 10.5 GeV C.M.E. 


D+ K-TT+n 


f 




D+X) = 640±14±35±43 


D° ^ K-TT+ 




a{e^e~ - 


D^X) = 1,521±16±62±36 






cj(e+e^ - 


D^X) = 1,579±55±102±63 


D*+ D^TT+ - 


(if-7r+)7r+ 


cr(e+e~ 


D*+X) = 583±8±33±14 


D*+ D^TT+ - 


(i<r^7r+7r+7r~)7r+ 


cj(e+e^ - 


D*+X) = 572±26±45±24 






cj(e+e^ 


D*^X) = 559±24±35±29 


D*o ^ 1)0^0 ^ 


(K-7r+7r+7r-)7rO 




D*^X) = 616±32±62±39 



We calculate < Xp > for the spectrum and for the spectra of D", D*^ and D*^ 
averaged over the decay modes. We supplement the data spectrum in the lowest bins using 
the JETSET spectra normalized to the spectra. We take the errors on these "borrowed" 
cross sections to be roughly comparable to the data in nearby bins. The results are shown 
in Table US 
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TABLE IV: < Xp> for the four charm mesons considered. The first error is statistical, the second 
systematic. 



Meson 


< Xp> 


Meson 


< Xp > 


D+ 


0.582 ± 0.008 ± 0.004 
0.570 ± 0.005 ± 0.004 


D*+ 


0.611 ±0.007 ±0.004 
0.596 ± 0.009 ± 0.004 



IX. OPTIMIZATION OF JETSET PARAMETERS 

Largely for internal use of our collaboration, we perform a simple fit of the spectrum 
(from the D^^K~7t^ decay mode) varying the three JETSET parameters that are most 
important for the shape of the spectrum. The first and second are the parameters a and b 
appearing in the "Lund Symmetric Fragmentation Function" j^, 2^ |: 



(1 — zY 
f{z) = N^ ^exp 



-b ■ 



(3) 



where z is the reduced energy xe, or momentum Xp, of the hadron and = +p'±, with 
m being the hadron mass and p± the component of the hadron momentum perpendicular 
to the jet axis. 

The third parameter is the probability Py that a meson of given flavor be generated as a 
vector meson, rather than pseudoscalar or tensor, Py = V/{P + V + T). The data indicate, 
as expected, that the majority of D°'s are not produced directly in the fragmentation of the 
charm quark, but from the decay of -D*'s. In JETSET these parameters are PARJ(41), 
PARJ(42) and PARJ(13). 

The result of the fit of the spectrum (in the K~'ir~^ decay mode) is: 

a = 0.178 ±0.007, 6 = 0.393±0.006, = 0.627 ± 0.015. 

Keeping Py fixed at the naive value Py = 0.75, we obtain a = 0.223±0.009 and b = 
0.438±0.005. In both cases the quoted errors are simple statistical errors. Correlation 
between parameters are not evaluated. The spectra resulting from these parameterizations 
are shown in Fig. d d HH |221 

Notice that we do not consider our results of D^, D*~^ and D*^ spectra in the optimization 
process. However, a posteriori we see, visually from the figures, that the spectra generated 
with these parameters seem to reproduce rather accurately also the D^, D*^ and D*^ ex- 
perimental distributions. However, it is not obvious which one of the two sets, the one with 
Py = 0.672 or the one with Py = 0.75, should be preferred. Furthermore, these parameters, 
while useful for the Monte Carlo simulation of D and D* spectra at the cm. energy of 
our and similar experiments, should not be taken as having general validity and theoretical 
significance. In fact, the D'^ spectrum generated by JETSET with our fitted parameters 
disagrees appreciably with the spectrum measured by the CLEO [3J| and BaBar 3^| col- 



laborations. It should be noted that the effect of these parameters may also be influenced 
by the value of other JETSET parameters. 



27 



X. CONCLUSIONS 



We have measured the momentum distribution of D^, D^, D*^ and D*^ produced in non- 
resonant e^e~ annihilation at a CME of about 10.5 GeV. These distributions can be used 
to guide and check QCD calculations of fragmentation functions needed to predict heavy 
meson production in both e~^e~ annihilation and hadron collisions at very high energy. The 
D° spectrum was used to determine the JETSET parameters that best reproduce it, and we 
found that, with these parameters, the D*, D*~^ and D*^ spectra (but not the spectrum) 
are also well reproduced. 
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APPENDIX A: PLOTS OF DETECTION EFFICIENCIES VS Xp 

In the following figures we show the detection efficiency dependence on Xp for all the 
mesons and decay modes analyzed. The detection efficiencies obtained from the signal and 
generic MC simulations are plotted, together with the curve resulting from the fit of their 
weighted average to a polynomial. 
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FIG. 23: Direct comparison of the detection efficiencies from signal and generic Monte Carlo and 
the result of smoothing their average: (a) for the D^^K^tt^tt^ channel, (b) for the D^^K^-k^ 
channel, and (c) for the ti~^ ti^ tt^ channel. 
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FIG. 24: Comparison of the detection efficiencies obtained from tlie signal and generic Monte 
Carlo and their smoothed average: (a) for the D*~^ — > D^7r~^ — > (i^~7r"'")7r+ channel, (b) for the 
D*+ ^ £)0^+ ^ (^-^+^-^+)^+ channel. 
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APPENDIX B: TABLES OF DIFFERENTIAL CROSS SECTIONS 

In the following tables, we report the quantity da/dxp in pb. Notice that the systematic 
and total errors are errors on the bin content {i.e., the first column). The first column 
of systematic errors is obtained from the rms spread of yields for the different procedures 
used to calculate the spectrum. The second column of systematic errors is derived from the 
"generic MC check" described in Sec. IVI Al These are the errors relevant to the shape of 
the spectra, i.e., they do not include the systematic errors that are common to the whole 
momentum range and that contribute to the error on the cross section (Sec. IV111|) . 
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TABLE V: da{e+e- D+X)/dx„ in pb: (D+^K-7r+TT+) 





/V/T / nif 




Errors (pb) 




Xp 






Systematic 


Tntnl 


U. iO-U.ZU 


iDi 


/ 


27 


3 


00 


n 9n n 91^ 

u.zu-u.zo 




i U 


53 


5 


Q9 


n 9^ n '^n 




70 


59 


6 


Of) 


vj.ov~yj.oo 


^lo 




68 


7 




n '^^ n zLn 

U.O0-U.4tU 


oyo 


00 


11 


11 


fin 


n zLn n zL^ 

U.4tU-U.4tO 




^9 
OZ 


14 


15 


OD 


U.4tO-U.OU 


1 nA9 

lU^Z 




16 


17 


00 




1 971 




20 


21 


00 


n n fin 


iOO i 


99 


21 


22 


00 


0.60-0.65 


1370 


19 


21 


22 


36 


0.65-0.70 


1291 


17 


20 


21 


34 


0.70-0.75 


1129 


15 


17 


18 


29 


0.75-0.80 


952 


13 


15 


16 


25 


0.80-0.85 


694 


10 


11 


11 


19 


0.85-0.90 


449 


8 


7 


7 


13 


0.90-0.95 


223 


5 


3 


4 


7 


0.95-1.00 


74 


3 


1 


1 


4 



31 



TABLE VI: daie+e- D^X)/dxp in pb; {D^^K-tt+). 
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TABLE VII: da{e+e- D^X)/dxp in pb; (£»°^i^-7r+7r+7r-). 
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TABLE VIII: da(e+e- ^ D*+X)/dxp in pb; (D*+^(is:-7r+)7r+). 
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TABLE IX: da{e+e- D*+X)/dx.„ in pb; {D*+^{K--k+tt+tt-)tt+). 
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41 

4i:± 


29 


22 


00 


60-0 65 


1276 


36 


30 


23 


52 


0.65-0.70 


1221 


32 


29 


22 


49 


0.70-0.75 


1096 


29 


26 


20 


44 


0.75-0.80 


915 


25 


22 


17 


37 


0.80-0.85 


715 


21 


17 


13 


30 


0.85-0.90 


533 


18 


13 


10 


24 


0.90-0.95 


317 


14 


8 


6 


17 


0.95-1.00 


122 


10 


3 


2 


11 



TABLE X: da{e^ e' D*^X)/dxp in pb; (L»*°^(iC-7r+)7r°). 





da/dxp 




Errors (pb) 




Xp 


(pb) 


Statistical 


Systematic 


Total 


0.15-0.20 


108 


121 


6 


2 


121 


0.20-0.25 


290 


121 


17 


7 


123 


0.25-0.30 


376 


112 


23 


9 


114 


0.30-0.35 


425 


104 


26 


10 


107 


0.35-0.40 


580 


95 


35 


13 


102 


0.40-0.45 


601 


92 


36 


14 


100 


0.45-0.50 


946 


99 


57 


22 


116 


0.50-0.55 


1061 


69 


30 


24 


79 


0.55-0.60 


1124 


61 


31 


26 


73 


0.60-0.65 


1186 


60 


33 


27 


73 


0.65-0.70 


1125 


56 


31 


26 


69 


0.70-0.75 


992 


48 


28 


23 


60 


0.75-0.80 


822 


47 


23 


19 


55 


0.80-0.85 


662 


36 


18 


15 


43 


0.85-0.90 


425 


28 


12 


10 


32 


0.90-0.95 


271 


24 


8 


6 


26 


0.95-1.00 


107 


22 


3 


2 


22 



34 



TABLE XI: da{e+e- D*°X)/dxp in pb; (L»*°^(K-7r+7r+7r-)7r°). 









Errors (pb) 




Xp 






Systematic 
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71 


28 


1 /LI 


0.60-0.65 


1367 


100 


75 


29 


128 


0.65-0.70 


1418 


85 


78 


30 


119 


0.70-0.75 


1235 


70 


68 


27 


101 


0.75-0.80 


954 


56 


52 


21 


80 


0.80-0.85 


764 


46 


42 


16 


64 


0.85-0.90 


581 


36 


32 


12 


50 


0.90-0.95 


317 


26 


17 


7 


32 


0.95-1.00 


131 


18 


7 


3 


20 



35 



TABLE XII: Differential cross sections da/dxp in pb for , D^, D*^ and D*^. The last three 
columns are weighted averaged over the two decay modes. The errors are the quadratic combination 
of the statistical and systematic errors, excluding the errors, discussed in Sec. IVIIII that affect the 
total cross section but not the shape of the spectrum. 



Xp 


D+ 











10- 


-0 


15 




i ( •) _l — L\jU 









15- 


-0 


20 


161 ±83 


431 ±186 




108 ±121 





20- 


-0 


25 


QOf) +Q2 


Koq +onQ 


J-'iU —LOW 


oqo +1 1 5 

Aj^J I 11(1 





25- 


-0 


30 
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882 ±84 


253 ±60 

^ijtj —i—\J\J 


387 ±111 
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-0 


35 


A-\ o _|_q4 


1 1 100 

±X(J1J _i_oo 


Q40 ICQ 


'-X.LiXj _l VXJKJ 
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-0 


40 


«qQ ICQ 


iRvn _|_q4 




c;q4 _|_qc 
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40- 


-0 


45 


qnq +5R 


Z. vj r .7 _l L iL/ 


624 ±46 







45- 


-0 


50 


1042 ±53 


2822 ±122 


920 ±50 


897 ±108 





50- 


-0 


55 


1271 ±38 


3194 ±56 


1108 ±32 


1085 ±70 





55- 


-0 


60 


1357 ±38 


3475 ±58 


1244 ±33 


1162 ±65 





60- 


-0 


65 


1370 ±36 


3371 ±56 


1286 ±32 


1230 ±64 
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-0 


70 


1291 ±34 


3007 ±51 


1248 ±31 


1198 ±60 
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-0 


75 


1129 ±29 


2549 ±46 


1113 ±29 


1055 ±52 
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-0 


80 


952 ±25 


2008 ±39 
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865 ±45 
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-0 


85 


694 ±19 


1383 ±30 


723 ±21 


694 ±36 
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-0 


90 


449 ±13 


829 ±21 


531 ±17 


471 ±27 





90- 


-0 


95 


223 ±7 


339 ±11 


310 ±12 


289 ±20 





95- 


-1 


00 


74 ±4 


90 ±5 


119 ±7 


121 ±15 
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